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Abstract 

The Grover search algorithm is one of the two key algorithms in the field of quantum computing, 
and hence it is of significant interest to describe it in the most efficient mathematical formalism. We 
show firstly, that Clifford's formalism of geometric algebra, provides a significantly more efficient 
representation than the conventional Bra-ket notation, and secondly, that the basis defined by the 
states of maximum and minimum weight in the Grover search space, allows a simple visualization 
of the Grover search as the precession of a spin-i particle. Using this formalism we efficiently solve 
the exact search problem, as well as easily representing more general search situations. 
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I. INTRODUCTION 

The Grover search algorithm jl-3], seeks to evolve a wave function, from some starting 
state \a), into the solution state \m), representing the set of all solutions, which upon 
measurement, will yield one element from this set j^-lg]. In order to analyze this evolution, 
typically an orthonormal basis \m) and Im -1 ) is defined, as shown on Fig (0Q) upon which 
the starting state \a) is plotted. However in this paper we use an alternative basis, the 
states of maximum and minimum weight, which allows the initial sate \a) and the solution 
state \m) to be symmetrically positioned in this space, allowing the conceptualizing of the 
Grover search process, analogous to the precession of a spin-| particle in a magnetic field, 
orecessing in this case, from the direction of the initial state \a) to the solution state \m) 
7|. This approach is similar to an SO(3) picture which has previously been developed jg], 
which also plots the path of the state vector during the application of the Grover operator. 



Clifford algebra has also been applied previously to Grover's algorithm |9Nl lj . however the 
approach adopted here combines the benefits of an efficient representation as well as an 
integral geometric visualization. 

A. The standard Grover search 

Given a search space of N elements, with M of these elements being solutions to a search 
query as defined by an oracle f(x), where by definition, f(x) — 1 if x is a solution, and 
f(x) = if x is not a solution, we set up two equivalence classes defined by the following 
two states: 



M ,, 

allowing us to define a uniform superposition starting state in terms of these two states as 

1 . IN — M [M, . 

i*> = 7/^ E w = V ^\H m > + v n |m) - (2) 

Grover's solution to the search process 4j, involves iteratively applying a unitary operator 
G defined by 

G = -G a G m = -(/ - 2\a)(a\)(I - 2\m)(m\). (3) 
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FIG. 1: Geometry of starting state \a). 

This operator applied to the n = [logiV] qubits representing the search space, rotates the 
state vector an angle given by 



sin 



e m 
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at each application, and after 



R< 



7T N 
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iterations, the wave function will lie close to the solution state \m). 



B. Modified basis vectors for the search space 



The Grover search space, is found to be isomorphic to an su(2) space [121 ]. and using 
the states of maximum and minimum weight, we can find a geometric picture of the search 
process, in a real three-dimensional space. Working from the two-dimensional complex space 
shown in Fig (pQ), using the well defined states \m) and |cr), we have four possible operators: 
\m)(a\, \a)(m\, \m)(m\ and |cr)(cr|. From these we define: 
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where a = (a\m), which may be complex, (3 = a/1 — |«| 2 , and i = \f—\ with 

(\a}(a\ - \m)(m\f 

IP 

We then find 



(9) 



[Ji,Jj] = ie ijk J k (10) 
P 

[Ji, Jj]+ = 3ij~2 

[K, Ji] = 0, 

where 5 is the Kronecker delta symbol, and e is the Levi-Civita symbol, confirming we have 
an su(2) algebra. Squaring the generators we find 

4 = Ji = Ji = ^(HH - HH) 2 = \p. (ii) 

We can easily check, P\o) = \cr) and P\m) = \m), with P 2 = P and the Casimir invariant 

c = jI + jI + jI = \p, 

which corresponds to a spin ~ system. We have the raising and lowering operators 

J± = Ji± i J 2 

and requiring J + \ f) = and J_| 4-) = 0, we find the states of highest and lowest weight 

9(9 9 \ 

It) = sec- \sm~\m) - e i5 cos J (12) 

|4.) = sec - I cos -|m) — e sm-|cr) 

where sin | = |a| and a = |a|e l5 . We then find J$\ t) = +|| t) an d ^3 1 I) — — ^1 4^}> as 
expected for a spin— | system. Writing \a) and |m) in this new basis we obtain 

\a) = e- i5 ^-cos^| t)+sin^U)J (13) 

, , f- cos - 
4' " 4 1 



■ , . 9 1 . 1 . 

|m) = — sm -| j) + cos -| 4-/ 



m ) = cos - T) + sin - !)■ 
Using these results, we can substitute back into the Grover iteration Eq. (j3J) to find 

G = -J + 2cos 2 ^| t>(t I +sin#| t)(| | -sin0| |)(t | + 2 cos 2 || 4,)(| |, (14) 
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which can be immediately written in matrix form as 

cos 9 sin 9 

G 



(15) 



- sin 9 cos 9 

which shows as expected, that the Grover operation rotates the state vector by an angle 9, 
where the starting state will be for this basis 



C. Clifford's algebra of three-space 



cos i 



sin 



(16) 



Using the orthonormal basis | f), | \), defined in Eq. (|T2|) . we now model the search 
process using the real associative algebra of GA. We define unit algebraic elements e±, e 2 , e 3 , 



such that e\ = e\ 



1, and for distinct i and j we have the anticommutation rule 



—efti [13(. The algebraic elements ei,e2,e3, define a three-dimensional space, and so 



we can define two vectors a = a\t\ + 0262 + 0363 and b = b\e\ + b 2 e 2 + 6363, then using the 
distributive law of multiplication over addition we find 

ab = (aiei + a 2 e 2 + a 3 e 3 )(&iei + b 2 e 2 + b 3 e 3 ) (17) 
= ai&i + a 2 b 2 + a 3 b 3 + (a 2 b 3 - b 2 a 3 )e 2 e 3 + (ai& 3 - a 3 6i)eie 3 + (a\b 2 - b\a 2 )e x e 2 
= a-b + taxb = a- b + aAb, 

where a - b is therefore the conventional dot or inner product and a A b is the wedge or outer 
product. In three dimensions we have the relationship with the conventional vector product 
that a A b = — ta x b, where we have defined the trivector 1 = eie 2 e 3 , which represents a 
signed unit volume. 

Using the product defined Eq. (fTTj) . with orthonormal basis elements, we find 

(18) 



C{Cj C{.Cj ~\~ Ci A Cj 6{j ~\~ ^ijk^k 



indicating that we have an isomorphism between the basis vectors ei,e 2 ,e 3 and the Pauli 



matrices through the use of the geometric product. We find that l 2 = t\t 2 e 3 e\e 2 e 3 



-1 



and we also find that 1 commutes with all other elements of the algebra and so can be used 
in place of the unit imaginary i = y/—T. The bivectors also square to negative one, that is 
[titj) 2 = ( e i e j)( e i e j) = ~ e i e j e j e % — — 1) which are used to define rotations in the plane of 
the bivector. 
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1. Rotations in 3-space with geometric algebra 

The Grover search process involves the incremental rotation of the state vector and, in 
geometric algebra, in order to rotate a vector by an angle |a| about an axis given by the 
vector a we define a rotor 

R = e ~' a/2 = cos(|a|/2) - i^- sin(|a|/2), (19) 

which acts by conjugation to rotate vector v = V\e\ + v^i + "^263, using 

v' = RvR) = e- ta/2 ve ta/2 . (20) 

The f is also called the reversion operation, which flips the order of the terms and switches 
the sign of 1. The bilinear transformation used to calculate rotations is preferred, because 
it applies completely generally, to not only rotating vectors, but also any components of the 
algebra, and also in any number of dimensions. 



2. Representing quantum states in Geometric Algebra( GA ) 



We can identify a simple 1:1 mapping from complex spinors to the bivectors of GA as 
follows 



13Hl5| 



zi\ t) + z 2 \ I) 



a + ia 3 
— a 2 + iai 



<B- ip = a + a k ie k . 



(21) 



This maps spinors to the even sub subalgebra of GA, which is closed under multiplication. 
Converting the complex spinors defined in Eq. ( Tl3j) . we find using Eq. ( 12T1) 



0.0 

\a) 1 — y —cos sm-te 2 = 

1 ' 4 4 


\m) 1 — y — sm - — cos -te 2 = 
4 4 

1 U 0.0 

\m) !->■ cos sm-te2 = e 

4 4 



3 ie 2 (5r/2-0/4) 



(22) 



-Le 2 e 



-Le 2 6/4 



We can now transform GA type spinors into a real space representation through the trans- 
formation 

S = (23) 
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which gives us the three-space vectors 



a = e ^/% e ~^ = e «»'/* ei = _ sin | 6l + cos ^ 3 



(24) 



m = e tea(»r/2-0/4) -tea(»r/2-0/4) 



-e^ e ' 2 = e^- 9 ^e 3 = - sin Q -e x - cos ^e 3 



m 



Hence the vectors a, m and m -1 , can now be plotted in a real Cartesian space analogous to 
the Bloch sphere representation as shown in Fig (j2J), and can be compared with [8]. As can 
be seen we use the e 3 (z-axis) from which to measure the angle 9, and is measured from 
e 1 . 



e 3 m Gm(T 




ei 



rn 



FIG. 2: Grover search in three-space based on states of maximum and minimum weight. The two 
possible precession axes are now e\ and &2i in order to rotate a onto m. 



II. THE GROVER SEARCH OPERATOR IN GA 

The action of the oracle G m on \m) is (/ — 2\m) (m\)\m) = — \m), which is to flip the 'm' 

n 

coordinate about the (m- 1 ) axis [5(. Reflections are easily handled in GA, through double 
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sided multiplication of the vector representing the axis of reflection, the action of the oracle 
being therefore 

m ± am ± = mam. (25) 

Using the starting state defined in Eq. ([21]) we find the action of the oracle on the starting 
state a as 

mam = ^9/2^-^9/^-^9/ 2^ = ^30/2^ = cog & + ^ 30 (2g) 

which is the required vector; c.f. Fig (j2J). 

The action of the other half of the Grover operator G a = I — 2\a)(a\ also produces a 
reflection, but this time about the a vector, which therefore implies a full Grover iteration 
of the starting state will be a(mam)a = (am) a (ma) = GaG^ , using associativity, giving 
the combined Grover operator as 

G = -am = e Le2e/2 e 3 e 3 e Le * e/2 = e Le2 °, (27) 

which is a significantly more compact form for the standard Grover operator, in comparison 
to Eq. Q. We can see by inspection, that the Grover operator represents a rotation of 26 
about the e2 axis, which will clearly rotate the vector a onto m, after an appropriate number 
of operations, as shown in Fig. [21 Hence, after k iterations we require the Grover operator 
G to rotate the vector a, defined in Eq. (1241) . onto the solution vector m, and so we require 

G k v& k = e ^V e2e / 2 e 3 e- tfce29 = e Le ^ 2ke+e ^e 3 = m, (28) 

using (e Le2d ) k = e Lke2e , and with m defined in Eq. ( |24l) . we therefore require 

^9+9/2)^ = ^(n-9/2)^ (2Q) 

and so by equating exponents, ignoring rotations modulo 2ti, we find the condition 

2k e+ e - = ir- e ~, (30) 



or 



and using 6 = 2 arcsin J ^ for a database with M solutions, we find 



' 71 lN (32) 



4 arcsin Ji 2 lV1/ 



the well known result for the standard Grover search. However, clearly, there is no guarantee 
that the formula will return k as an integer, and because it will need to be rounded to the 
nearest integer describing the number of Grover operations, then we will not always return 
exactly the solution space upon measurement. However we can modify the search slightly, 
in order to guarantee that k will be an integer, and hence reliably return the solution state 
\m). 



A. Exact Grover search 

The Grover operator, defined in Eq. ([3]), can be modified so that it rotates the starting 



state \a) exactly onto the solution states |m), thus finding a solution with certainty 12l. Il6|. 
In order to create the exact Grover search, the Grover operator is typically generalized to 

G = _(/ _ (i _ e^)\a)(a\)(I - (1 - d+>)\m){m\), (33) 



so that when the oracle identifies a so 
function and not just the scalar —1 



ution it applies a complex phase e 1 ^ 2 to the wave 
12] . This has the effect of slightly slowing down the 
search process, but it then allows the solution state \m) to be reached exactly using an 
integral number of iterations. 

A reflection can be viewed as a rotation by tt in one higher dimension, so if we rotate by 
an angle 02 about the m axis, which will be clockwise as viewed from above the e% axis, we 
obtain the oracle 

Q m _ e ^(sin(0/2)ei+cos(0/2)e 3 )_ (34) 

For (f>2 — tt we find G m = i(sin(#/2)ei + cos(9/2)e 3 ) = im, so that the action of the oracle 

GmoG^ = Lma(—Lm) = mam, (35) 
which gives the same result as the standard Grover oracle found previously Eq. f l25|) . similarly 

Gfj = e -^(-sin(e/2)e 1+ cos(0/2)e3) (gg) 

will be a rotation about the a axis. Hence for the exact search for the Grover operator we 
have 

G = —G a G m (37) 

_ _ e -^(-sin(e/2)e 1 +cos(e/2)e3) e ^(sin(e/2)ei+cos(e/2)e 3 ) 



which when expanded gives 



G = cos — cos h cos U sm — sin h sm sin -tei (do) 

2 2 2 2 2 2 y ' 

. 01 . 02 . Q 6* . 01 - 02 

+ sm — sm — sm ULe 2 — cos - sm te 3 , 

2 2 2 2 

for the general Grover operator, which is now interpreted as a rotation about a general 

precession axis. 



1. Phase matching 

We can see from Fig. |2j which uses the alternate orthonormal basis | f) and | 1), that a 
and m now lie in the plane of e\ and e^, and hence using a geometric argument the Grover 
precession axis must therefore lie in the plane of e\ and e 2 in order for the a vector to be able 
to be rotated precisely onto the m vector. Hence we need to eliminate the component in 



the precession axis, and hence, by inspection of Eq. ( IdTj) . we require 0i = 02, which is the 
well known phase matching condition [lfl [l^] . Hence the exact search will be in the form 

Q _ e t/3(sinaei+cosoe2) (39) 

where we find 

sm — = sm - sm — (40) 

2 2 2 V ; 

6 
cot a = cos -tan—, 
2 2' 

which can be re-expressed assuming a normalization factor Z as 

Q _ gt/3(cos |ei+cos I sin ^e 2 )/Z {^) 

which shows clearly the precession plane perpendicular to the vector cos te% + cos | sin 
and if we select = 7r, we recover the standard Grover search operation. 

To calculate for the exact search we, once again have the vector equation, given by 
Eq. (I2"g|) . and substituting our modified Grover operator, along with Eq. (|2"4"|) . we find 

e ifc/3(sinaei+cosae2) e te 2 e/2 e ^ e — tfc/3(sin«ei+cosae 2 ) _ g — Le 2 0/2^ (42) 

which can be rearranged to 

gi./c/3(sinaei+cosae2)gte 2 9/2g;,fc/3(sinaei+cosoe2)gte20/2 _ ^ (43) 
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or 

Now, because we can always replace two consecutive precessions, with a single precession 
operation, we can write 

for some unit vector v. Thus, from Eq. (144 p . we need to solve 

(e"*) a = e 2tK " = cos 2k + 6v sin 2k = - 1 (46) 

and so clearly k = |. Thus the right hand side of Eq. f j4"5l) . is equal to tv, implying that the 
scalar part is zero. Expanding the L.H.S. of Eq. ( I4"5l) . and setting the scalar part to zero, 
we find 

6 9 

((cos k{3 + l sin kf3(sin ae\ + cos ae-^j) (cos - + i sin -e2))o (47) 

= cos kB cos sin kB sin - cos a = 0. 

y 2 2 

Re-arranging this equation we find 

sin ^ 

cot k/3 = tan - cos a = — (48) 

2 ^cos 2 | + cot 2 | 

Isolating fc, we find 

arccot 



z arcsm(sm | sm 

Using calculus we can find the minimum at = ir, which thus returns the number of itera- 
tions for the standard Grover search given by Eq. (I3TT) . which shows that this modification 
fails to speed up the search 18N22|. However we are able now to set <p in Eq. fT49|) . so as 
to make k an integer, which will therefore be the fastest exact search possible. Hence the 
minimum integer iterations will be 



k,, 



7T 1 

26 ~ 2 



(50) 



Substituting back into Eq. (I49p and re-arranging we then find an expression for 



, 9 d>\ I sinf , 

2fc m arcsin I sin - sin — ] = arccot — , (51) 



2 2 



'cos 2 | + cot 2 2 
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which we can simplify to give explicitly 

. 6 . 7T 9 

sm — = sm — esc — . (52) 

2 4:k m + 2 2 v ; 

We have <fi now determined directly from the known 9 and k m defined in Eq. (j3J) and Eq. fl50|) 
respectively, thus solving the exact search using the Grover operator defined in Eq. ( I39lh 

An example using this formula for the exact search is given in the Appendix, which shows 
how the starting polarization vector now rotates exactly onto the solution states, as required. 



B. General Exact Grover search 

Most generally we can write the Grover operator as 

G = -UI~,U- l G m (53) 

where 

/ 7 = J + (e^-l)|7>< 7 | (54) 
where we normally choose 7 = |0) = |0 . . . 0), 0, 23, 24]. For U = H we have 



G = -UI 1 U- l G m = -(/ + (e^ 1 - 1)#| 7 )( 7 |# = -(/ + (e^ 1 - l)\a)(a\ = -G a G m . (55) 

So with this modified operator we effectively use a modified vector to a, namely the vector 
7 = U\0), giving 

| 7 ) = -e-^cos^l t)+e^ /2 sin^||), 

equivalent to a starting polarization vector 

• ^° a ■ @° ■ a 1 ^° (za\ 

7 = — sm — cos 0oei — sm — sm o e2 + cos ~^ e 3- (ooj 

Comparing this with the polarization vector for the standard Grover search a = — sin |ei + 
cos |e3, as shown on Fig [2], we see that we have changed the projection in the e% direction 
by changing 9 to 9q, and hence rotated the vector in the plane given by the vector <p . 
If 0o = 0, then we recover the standard exact Grover search. As this is a unit vector, we 
simply adapt G G to rotate about this new vector, that is we have 

Q _ g— (.7^1/2 _ g— sin -Q cos^oei-sin ^ sin0 o e2+cos ^-ea) (57) 
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and hence for the general exact search we have 




ni 



— sin cos 4>oei— sin sin </>oe2+cos -^-ez) i</>2/2(sin(0/2)ei+cos(0/2)e3) 



(58) 



However, as a more elegant alternative, we can simply adjust our basis states, given by 
Eq. ffl2|) . and then the exact solution, given by Eq. f )39|) . immediately applies. 

III. SUMMARY 

The Grover search algorithm is a central algorithm in the field of quantum computing, 
and hence it is important to represent it in the most efficient formalism possible. The two 
main strengths of geometric algebra are its method of handling rotations and its integral 
geometric representation, and hence its perfect suitability in describing the Grover search. 
We find Clifford's geometric algebra, provides a simplified representation for the Grover 
operator Eq. ( )27|) and a clear geometric picture of the search process. Using the states of 
maximum and minimum weight, we find that we can interpret the search process as the 
precession of a spin-| particle, thus providing a simple visual picture, as shown in Fig |2j 
This is not possible with the standard formalism as it requires two complex axes, forming a 
four-dimensional space, and hence difficult to visualize. We also find that the exact Grover 
search Eq. ( |39|) has an efficient algebraic solution, as shown in Eq. (152]) . Improved intuition 
obtained via the use of Clifford's geometric algebra, may possibly enhance the search for 
new quantum algorithms. 

Appendix A: Example of an exact search over 16 elements 

After k iterations we have the polarization vector 



p = G k aG ]k 



(Al) 



tfc/3(sin oei+cos ae2 „ te 20 /2 — tfc/3(sin«ei+cos ae2 



+ I cos - cos 2/3k — cos a sin - sin 2/3k e^. 




sin 2 a sin — h sin - cos 2 a cos 2(3k + cos a cos - sin 2(3 k I e± 
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For 16 elements we find from Eq. (!50l) k m = 3, and we then find <fi from Eq. (|52j) for an exact 
search of (f) = 2.19506. This gives the polarization vector after k iterations 

P = -(0.0546434 + 0.195357 cos 2/3A; + 0.855913 sin 2(3k) ei (A2) 
+ (-0. 10332 + 0. 10332 cos 2/3k + 0.452673 sin 2(3k)e 2 
+ (0.968246 cos 2/% - 0.220996 sin2/3fc)e 3 . 

Using a and defined in Eq. (|4"01 . beginning from a starting vector a = (—0.25,0,0.9682), 
with a required solution vector m = (—0.25, 0, —0.9682), we generate a polarization vector P 
as 

a = (-0.25,0,0.9682) (A3) 
Ga = (-0.8456,0.315,0.4309) 
G 2 a = (-0.8456,0.315,-0.4309) 
G 3 a = (-0.25,0,-0.9682) 

thus producing the exact solution m after k m = 3 iterations as required. 
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